Rydberg Borromean Trimers 



00 

o 
o 

(N 



Ivan C. H. Liu. Jovica Stanojcvic. and Jan M. Rost 

Max Planck Institute for the Physics of Complex Systems, 

Nothnitzer Strasse 38, D-01187 Dresden, Germany 

(Dated: July 13, 2008, version 1.6) 

A Rydberg and a ground-state atom can form ultralong range diatomic molecules provided the 
interaction between the ground-state atom and the Rydberg electron is attractive [C. H. Greene, 
et al., Phys. Rev. Lett. 85, 2458 (2000)]. A repulsive interaction does not support bound states. 
However, as we will show, adding a second ground-state atom, a bound triatomic molecule becomes 
possible constituting a Borromean Rydberg system. 

PACS numbers: 36.40.Gk, 31.15.Qg, 36.40.Wa 
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With the advent of ultracold atomic physics new forms 
of atomic systems become possible. The trilobite, ar- 
guably one of the most exotic molecules, is a ground- 
state atom bound to a Rydberg atom by the polarization 
potential of the latter. The existence of this extremely 
long range molecule was proposed, along with its name, 
by Greene et al. |l|. This idea generated considerable 
interest [3, S S H lal i although a clear experimental veri- 
fication has not been achieved so far but may be possible 
within a BEG making use of its density orders of mag- 
nitude larger compared to an ultracold gas [31 ■ More- 
over, previously unexplained satellite peaks in thermal 
spectra of molecules could be traced to the formation of 
such long molecules Q . The basis of the molecular bond 
in such systems is a negative scattering length leading 
to an attractive atom-Rydberg atom potential which is 
in addition modulated by the electronic Rydberg wave- 
function. An extension of the trilobite to polyatomic 
Rydberg molecules involving several ground-state atoms 
and one Rydberg atom was investigated in Q. There, it 
was found that certain geometric arrangements have adi- 
abatic molecular potentials which are many times deeper 
than in the diatomic case. 

However, a Rydberg molecule can even form with a 
positive scattering length, i.e., repulsive interaction, as 
we will show in the following. While a diatomic molecule 
is indeed not possible under repulsive interaction, a tri- 
atomic molecule can be formed which we call Borromean. 
The name comes from a similar phenomenon in nuclear 
physics Q- The essence of a Borromean system is the 
existence of a three-body bound state although none of 
the two-body subsystems are bound. For the origin of 
"Borromean rings" see [i3| • 

In the following, we will show that the building princi- 
ples of such ultralong range molecules allow the formation 
of a Borromean trimer using neon as an example. This 
is a molecule formed out of three linearly aligned neon 
atoms where the one in the middle is Rydberg excited. 
The mutual interaction of two neon atoms is repulsive 
and none of the atomic pairs has a bound state. Nev- 
ertheless, the triatomic system can have a bound state. 
We will use atomic units unless stated otherwise. 

Since the Rydberg wavefunction does not vary over the 
extension of a ground-state electronic wavefunction, the 



interaction can be described by a Fermi-pseudopotential 
lull, which amounts to a contact interaction 



yF(r;ri) = Li5{y-vi). 



(1) 



where r is the coordinate of the Rydberg electron and ri 
is the position of the ground-state atom A measured from 
the Rydberg atom. Furthermore, Li = — 27rtan(5o/A:(ri) 
is the energy- dependent scattering length of the collision 
system e~ + A containing the s-wave phase shift 5q and 
the local momentum fc(ri) of the Rydberg electron, re- 
lated to the distance ri through the energy conservation 
k{riY /2 = — l/2n^ -|- \/R. In the low-energy limit, the 
modified effective range theory [12] expresses Li analyti- 
cally in terms of the zero-energy scattering length Lq and 
the polarizability of the atom a. 
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For impact energy larger than 0.003 cV in the e~-|~Ne sys- 
tem, ah initio phase shift data |13| allow one to connect 
the scattering length of Eq. ^ to higher energies accu- 
rately, using the values Lq = 2.613 a.u. and a = 0.2218 
a.u. from [l3[. 

For N ground-state atoms at positions r j , the full elec- 
tronic Hamiltonian reads 
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where Vc{r) is the Goulomb interaction of the Rydberg 
electron with its mother ion. 

The eigenenergy of a dimer ft,2 may be obtained by di- 
agonalizing the interaction Vp within a degenerate mani- 
fold of hydrogenic eigenfunctions $„(m of principal quan- 
tum number n. This leads to the (unnormalizcd) solution 
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*^im(ri)*nim(r). 



(4) 
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Low-/ states may have different energies due to large 
quantum defects, so they need to be excluded from the 
basis formed with states of Goulomb energy — l/2n^. 
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with the normalization factors 
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FIG. 1: The adiabatic potential Eq. (|8} with ri — R for the 
(n = 30) neon Rydberg dimer (dash-dotteded black) and the 
symmetric stretch cuts Eq. (|14|) of trimer in solid blue (+) 
and dashed red (— ). 



Making use of the addition theorem for spherical har- 
monics {x is a unit vector) [lj|, we reduce Eq. (j4|) to 

V'„(r;ri)= ^ (t)niiri)(j>„iir)— — P/(rir), (5) 



/>;„ 



with the radial hydrogenic eigenfunctions (j^ni- It 
is convenient to specify the overlap intergral N12 = 
(-(/'n(r;ri)|'0„(r;r2))r which yields 

2/4-1 
^12 = ^0nz(r2)0n;(n)^^Pz(r2fi). (6) 

I 

Clearly, the normalized eigenfunction Eq. ^ reads 
vE'„(r;ri) = ifin{r;ri)/N^{ . For future reference we de- 
fine 



Tvfi = W'n(i-i; n) ± V'n(ri; "ri))/2 • 



(7) 



The relation Eq. ([6|) allows one to write the dimer energy 
E2{ri) = {^n\h2{ri)\'^n) + 1/2^2 in the compact form 



E2iri)^LiNn 



(8) 



The repulsive interaction Vp for the neon Rydberg dimer 
Ne-I-Ne* E2 is shown as a dash-dotted line in Fig. [1] 
Clearly, this system does not have bound states. 



However, we will show that adding another neon 
ground-state atom allows for bound states of the neon 
trimer with Hamiltonian 

2 
/i3 = y + Vcir) + Li5{v - n) + L2<5(r - V2) , (9) 

where Li > 0. The following basis functions are de- 
fined which are symmetry-adapted linear combinations 
of dimer eigenfunctions with gcradc/ungerade (±) sym- 
metry (since the two neon atoms 1 and 2 are identical), 

*±(r; n, r2) = [*„(r; n) ± *„(r; r2)]/A^±/' , (10) 



iV± = 2 ± 2 



A^i 
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V^11^^22 



(11) 



In this set of basis functions, the Hamiltonian h^, forms a 
2x2 matrix with diagonal terms E^^^ ~ (^±|/i3|'I'±) + 
l/2n^ of the simple form 



^d'ia(ri,r2) 



LiNn + L2N22 



N± . (12) 



The off-diagonal terms Eoff = (*±|/i3|*:f) + 1/271^ read 
LiNii + L2N22 



£'off(ri,r2) 



-0^^ 



(13) 



A symmetric cut through this potential for r2 = — ri = 
R is shown in Fig. [1] This cut has the appealing analyt- 
ical form 



^dl(ri, 
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(14) 



Hence, depending on gerade/ungerade symmetry only 
even/odd components I in the sum of Eq. ^ contribute. 
The symmetry-induced separation of basis functions with 
even and odd I results in a more pronounced undulation 
of the corresponding potentials (as can be seen in Fig.[T]), 
since there is less intereference from wavefunctions with 
different nodal structure. For the full solution the 2x2 
potential matrix can be diagonalized whose eigenvalues 
give the adiabatic potentials in the analytical form 



S"/^(ri,r2) = 



± 



LiNii + L2N22 



V(LiA^ii - ^2^22)2 + 4LiL2Ar22 



(15) 



To determine if the neon Rydberg trimer can support 
bound states we make a normal mode analysis on the 
electronic surface from Eq. pS]) . to which the outermost 
minimum belongs. To this end it is useful to define sym- 
metry adopted Jacobi coordinates for the nuclear motion. 



rA-(ri+r2)/2, R=(ri-r2)/2. 



(16) 



In these coordinates, the nuclear Hamiltonian for total 
angular momentum L = is given by 



H = 



1P2 



^Pa 
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E''^''irA + Rz,rA-Rz), (17) 



where m is the mass of a neon atom, and we have chosen 
the coordinates such that R = Rz along the body fixed 
z-axis. We may re-express the potential energy surface 
as^"/^ = f"/'^(i?,z,p). 

In the following, we quantize the trimer with n = 30 
in a separable approximation closely related to the nor- 
mal mode analysis for molecules. With the condition 
d£'^/dR = d£'^/dp = d£'^/dz = we determine the equi- 
librium points {R, p, z) = {Ri, 0,0) of the potential. Fig- 
ure [1] reveals many minima in the symmetric stretch (SS) 
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FIG. 2: The adiabatic potential Vs(p) (see Eq. Hill) at the 
outermost equihbrium point _Ro as a function of the bending 
coordinate p for the (n = 30) neon Rydberg trimer. The 
soUd blue and the dashed red correspond to E'^^^ and £"^1^^ in 
Eq. (|12|l respectively. The shaded area shows the quantized 
state. 
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FIG. 3: The adiabatic potential Va[z) (see Eq. dTS])) at the 
outermost equilibrium point 7i!o as a function of the asymmet- 
ric stretch coordinate z for the (n = 30) neon Rydberg trimer 
as well as the corresponding potential for the £" potential 
(upper curve). 



coordinate i?; we will quantize the outermost minimum 
at i?o = 850 au. As can be seen on the same figure, 
tunneling of the quantized states in the R coordinate de- 
pends sensitively on the form of the potential. Hence, 
instead of a formal harmonic normal mode expansion 
about the equilibrium point, we simply split the poten- 
tial in several contributions, leaving only one coordinate 
as a variable in each term while fixing the others at their 
equilibrium values. This allows for an explicit numer- 
ical solution in each degree of freedom under the total 
potential 

£\R, z,p) = Vs{R) + Va{z) + Vb{p) - 21/o , (18) 
with 

Vo = f''(i?o,0,0) (19) 

Va{z) = £^{Ro,0,z) 

Vs{R) = E+^{Rz,-Rz) 

Vb{p) = E;^.^{px + Roz,px- Roz). 

The modes correspond to molecular types of normal 
mode vibrations along the SS (in i?), the bending (in p) 
and the asymmetric stretch (AS) (in z) with the product 
wavef unction 



where 



^'>iR,z,p) = x{Rmz)ap)^ 



1 d^ 



AmdR^ 



(20) 



+ Vs{R) x{R) - EsxiR) (21) 



and similarly for the other coordinates. 



The quantized bending motion has a zero-point en- 
ergy at Eb = 1.04 X lO^^au as shown in Fig.[l The SS 
potential shown in Fig. [1] supports a ground state with 
energy Eg = 8.11 x 10~^ au and an excited state. Both 
of them are in principle resonances which decay through 
tunneling in the separable approximation. However, the 
ground state is practically stable. Last, we come to the 
AS which has near equilibrium (z = 0) the form of an 
inverted oscillator (Fig. [3]). The downwardly-bent curve 
falls on the left and right into deep but narrow wells, 
which arc the normal Ncj electronic ground-state po- 
tentials, to which the asymptotic polarization potential 
—a/2rf of the Rydberg excited level is smoothly con- 
nected. For even shorter distances, the ah initio data for 
the NeJ potential from Ref. [l^l are used. The potential 
landscape is similar to those of triatomic ABA molecules 
for which so called hyperspherical resonances exist [16| . 
The corresponding quasibound resonant eigenfunctions 
are highly excited with many nodes along the AS and the 
resonant energy Ea almost touches the top of the barrier, 
i.e., the equilibrium point Ea ~ Vq = 7.30 x lO^^au. To 
a good approximation the lifetime of the resonance can 
be estimated by the curvature kA of the potential at Vq, 
ta ~ 2'KUJ~^ = 2-KyJ ^A/kA = 57 ns with ^a = 2to/3. 
Intuitively, this curvature determines the time scale on 
which probability density "slides down" the inverted os- 
cillator potential towards the dimer diatomic decay chan- 
nel. To summarize, the total ground state energy of the 
Borromcan trimer measured from the Rydberg energy 
— l/2n^ is given by 



Ea = Es + EA + 2Eb - iVo = 1.43 x 10" 



au. 



(22) 



This means a blue shift of 940 MHz from the Rydberg 
excitation line and the lifetime is limited by the decay of 
the AS mode to about r ~ 57 ns. Due to the cylindri- 
cal symmetry of the collinear trimer in equilibrium, the 
bending motion in (p, 0) is doubly degenerate, hence it 
contributes 2Eb- 
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FIG. 4: The adiabatic potential for the lower surface of the 
(n — 30) neon Rydberg trimer in colhnear configuration. 
Symmetric (SS) and asymmetric stretch (AS) modes are in- 
dicated, as well as the dynamics of the trimer. 



For n = 30, the bending can be neglected due to the 
localized motion about p = and the long bond-length. 
Figure |4] gives an overview of the dynamics of the trimer 
in the plane of the SS and AS coordinate (in coUinear 
geometry with p — 0). The system rolls down from the 
top of the saddle point through a potential valley into 
the two-body region. 

The lifetime of the Rydberg Borromean trimer can be 
optimized by varying the Rydberg excitation n: Higher 
n leads to shallower potentials which enlarges the life- 
time in the AS but shortens it in the SS and the bending 
due to a stronger tunneling. Hence, there exists some n 
with the maximum r. To determine its value, we first ob- 
tained a general expression for the curvature of the saddle 
in the AS, which give s an estimate for the decaying time 
TA ~ n y^2TT-^fiA/LQ. This is obtained by approximat- 
ing the shape of the saddle with a parabola of height 
LQ/irn^ and width 2n^. We then numerically calculate 



the tunneling times ts and tb in the SS and the bending 
by solving the corresponding Schrodinger equations and 
matching the asymptotic phases, for increasing n at the 
outer-most minimum until the tunneling in the bending 
becomes too large. For Ne3, it was found that this opti- 
mal value occurs at about n = 55, where the overall life- 
time is as large as r = (I/T5 + 2/tb + l/r^)^^ = 0.46/iS. 

We note in passing that a similar analysis of the Ry- 
dberg trimer under negative scattering length for ®^Rb 
produces a bound state, red shifted by -12.7 GHz with 
respect to the Rydberg line, with the corresponding po- 
tential minimum Vq = —14.3 GHz at (865,0,0). 

One may question the validity of the adiabatic approx- 
imation in the presence of a series of conical intersections 
as obvious from the crossings in the diagonal part of the 
hamiltonian matrix Eq. (J14p and Fig. [TJ However, the 
vibrational ground-state wavefunction xo is strongly lo- 
calized about Rq, suppressing the effect of coupling even 
at the closest conical intersection. 

Future work will consider the effect of a series of conical 
sections in a situation where they are dynamically active, 
i.e., not sitting on a decaying background potential which 
renders the adiabatic wavefunction exponentially small 
under a larger barrier at the location of the intersections. 

To summarize, we have demonstrated that a Rydberg 
atom may form with two ground-state atoms an ultralong 
range Borromean molecule with finite but long lifetime 
for positive scattering length of the Rydberg electron- 
ground-state-atom interaction. We call it Borromean 
since the corresponding diatomic molecule Rydberg and 
ground-state atom, do not have a bound state or long 
lived resonance. 
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